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Abstract 

A search for nucleon decay into two-body hnal states containing K'^ 
mesons has been conducted using the 963 metric ton Soudan 2 iron track- 



ing calorimeter. The topologies, ionizations, and kinematics of contained 
events recorded in a 5.52 kiloton-year total exposure (4.41 kton-year fiducial 
volume exposure) are examined for compatibility with nucleon decays in an 
iron medium. For proton decay into the fully visible final states /U^Kg and 
e^Kg, zero and one event candidates are observed respectively. The lifetime 
lower limits {t/B) thus implied are 1.5 x 10^^ years and 1.2 x 10^^ years, re- 
spectively. Lifetime lower limits are also reported for proton decay into ^"^K^*, 
and for neutron decay into z^K^. 



Typeset using REVT^ 



I. INTRODUCTION 

A. SUSY predictions for nucleon decay 

Supersymmetry (SUSY) proposes the existence, presumably at high mass scales, of a new 
fermion or boson partner for each boson or fermion encompassed by the Standard Model. At 
the expense of doubling the number of elementary particles, SUSY Grand Unified Theories 
(GUTs) provide a natural solution to the hierarchy problem, allow the extrapolations of the 
running coupling constants (which have now been precisely determined at low energies and 
electroweak unification scales) to converge at a high energy scale, make accurate predictions 
of the Weinberg angle, and incorporate gravitation when locally gauged [0-^. 

Supersymmetric GUTs, such as SUSY SU(5), permit the same processes for nucleon 
decay as do their non-supersymmetric counterparts; however, the larger mass requirements 
for the GUT-scale lepto-quark bosons generally preclude the possibility of observing nucleon 
decay with lifetime divided by branching ratio r/B of less than ~ 10'^'^ years. Interestingly, 
other processes involving supersymmetric particle loops arise that can mediate nucleon decay. 
An example is shown in Fig. |l| wherein the decay is suppressed by the product of the color 
triplet Higgsino and Wino masses rather than by the supersymmetric version of the X and 
Y lepto-quark bosons PJ^. Nucleon decay diagrams of this type give integrals that vanish 
unless the transitions involve intergenerational mixing. Consequently, final states containing 
strange mesons are to be expected, with lifetimes predicted to be ~ 10'^^ — 10^^ years [|^§. 

Table I lists some proposed SUSY theories along with their predicted decay modes of 
leading branching ratios and their lifetime predictions. The expectation that two-body 
decays yielding a positive-strangeness K meson are predominant, is seen to be a common 
theme. 



B. Previous experimental searches for K^ modes 

Experimental searches for baryon instability in general and for nucleon decay into two- 
body modes containing K-mesons in particular, have been carried out for more than two 
decades. Among the early searches which included strangeness modes were the planar iron 
tracking calorimeter experiments of KGF |T^,n], NUSEX |]T3|, and Soudan 1 [W\. Lifetime 



lower limits and also observation of a few nucleon decay candidate events, were reported. 
Among candidate sightings was a three-prong event recorded by NUSEX having kinematics 
compatible with 

p^^+K°, K°^7r+7r-. (1) 

With more exposure however, this candidate was found to be consistent with background 
arising from interactions of atmospheric neutrinos within the NUSEX calorimeter. 

Lifetime lower limits for proton decay (|l|) (but including the K^ -^ 7r°7r° mode) and for 
neutron decay 

n^z/K°, K^^vrV. (2) 



were reported by the 1MB water Cherenkov experiment in 1984 |T^. These limits were 
improved upon and the scope of search extended to other K-meson modes in subsequent 
investigations by 1MB |lT8|-pO|. Lifetime limits were also obtained with the HPW water 



Cherenkov detector [^]. More recently, nucleon decay limits for many modes have been 
pubhshed, including the K° modes considered in the present work, based upon the 7.6 kton- 



year exposure recorded by IMB-3 22 



Nucleon decays involving K*^ and K"*" mesons were searched for using the iron tracking 
calorimeter experiment of the Frejus collaboration; lifetime lower limits are reported in 
Refs. [p^-p5|. Among the highest of published lifetime lower limits for lepton plus kaon 
two-body modes are the 1989 results by the Kamiokande water Cherenkov experiment 
based upon a 3.76 kton year exposure of KAM-I and KAM-IL 



A search for the proton decay mode p -^ v¥j^ in the Soudan 2 experiment has been 
reported in a previous pubhcation [^. A more stringent hfetime hmit for this mode has 
recently been pubhshed by Super-Kamiokande |^. In this work we report our search in 
Soudan 2 for nucleon decay into two-body modes involving K° mesons, namely proton decay 
into /i^K° and e+K° with K° — » K^ or K°, and neutron decay into z^K°. Soudan 2 has the 
capability to observe these leading SUSY decays if indeed their lifetimes are within the 
ranges indicated in Table I. 

II. DETECTOR AND EVENT SAMPLES 
A. Tracking calorimeter and active shield 

Soudan 2 is a 963 metric ton iron tracking calorimeter which is currently taking data. 
It is located at a depth of 2340 ft (2070 mwe) on the 27th level of the Soudan Underground 
Mine State Park in northern Minnesota. Data-taking commenced in April 1989 when the 
detector's total mass was 275 tons. The modular design allowed data-taking to proceed while 
additional 4.3-ton calorimeter modules were being installed. The detector was completed 
in November 1993 at 963 tons. The total (fiducial) exposure analyzed here is 5.52 (4.41) 
kiloton-years (kty), obtained from data-taking through October 1998. 

Each of the 224 modules that comprise the calorimeter consists of layered, corrugated 
iron sheets instrumented with drift tubes which are filled with an Argon-C02 gas mixture. 
Electrons liberated by incident ionizing particles drift to the ends of the tubes under the 
influence of an electric field sustained by a voltage gradient which is applied along the tubes. 
The drifted charge is collected by vertical anode wires held at a high voltage, while horizontal 
cathode pad strips register the image charges. Two of the position coordinates of the ionized 
track segments are determined by identifying the anode wire and cathode pad strips excited; 
the third coordinate is determined from the drift-time of the liberated charge once the event 



'start time" has been determined 29,30 . 



Surrounding the tracking calorimeter but mounted on the cavern walls and well separated 
from calorimeter surfaces, is the 1700 m^ active shield array comprised of 2 - 3 layers of 
proportional tubes [^ . The shield enables us to identify events which are not fully contained 
within the calorimeter. In particular, it provides tagging of background events initiated by 
cosmic ray muons, as will be described in sections [11 B| and |11D| below. 



In searches for nucleon decay into multiprong modes, the Soudan 2 iron tracking calorime- 
ter of honeycomb lattice geometry offers event imaging capabilities heretofore not achieved 
by the planar iron calorimeter experiments or by the water Cherenkov experiments. For the 
multiprong events of this work, Soudan's fine-grained tracking enables event vertex loca- 
tions to be established to within to 3 cm in the anode or cathode coordinate (drift-time 
coordinate) relative to the honeycomb stack (to the rest of the event). This spatial resolu- 
tion enables discrimination, based upon proximity to the event vertex, between a prompt e^ 
shower and a photon initiated shower. In Soudan 2, ionizing particles having non-relativistic 
as well as relativistic momenta are imaged with dE/dx sampling. Consequently, final state 
topologies are more completely ascertained, and proton tracks can be distinguished from 
charged pion and muon tracks. These capabilities are valuable for rejection of background. 

B. Event samples 

Events used in this analysis are required to be fully contained within a fiducial volume, 
defined to be the 770-metric-ton portion of the calorimeter which is everywhere more than 
20 cm from all outer surfaces. There are four distinct samples which are either simulated 
(Monte Carlo events) or collected (data events) : 

i) We use full detector nucleon decay Monte Carlo simulations, including detector 
noise, to generate samples of the various nucleon decay processes considered in this 
work. Hereafter we refer to these events as nucleon decay MC events. Details concern- 
ing generation and utilization of such events are given in Ref. |2^ . 



a) A full detector neutrino Monte Carlo (MC) simulation, which includes detector 
noise, is used to generate events representing all charged-current and neutral-current 
reactions which can be initiated by the flux of atmospheric neutrinos. These neutrino 
MC events are injected into the data stream prior to physicist scanning and are sub- 
sequently processed and evaluated using the same procedures as used with the data 
events. The Monte Carlo program is based upon the flux predictions of Barr, Gaisser, 
and Stanev [04|, ^ and is described in a previous Soudan 2 publication |^. Among 



the neutrino reactions simulated, care has been taken to include K /K channels; these 



O'TtO 



arise at low rate from AS = associated production and from AS = 1 processes ^6 
All of the atmospheric neutrino MC samples in this paper correspond to a flducial 
detector exposure of 24.0 kton- years. 

in) There are contained events for which the veto shield registered coincident, double- 
layer hits. Such events usually originate with inelastic interactions of cosmic ray muons 
within the cavern rock surrounding the detector. The muon interactions eject charged 
and neutral particles out of the rock and into the calorimeter; these particles sub- 
sequently interact to make events in the detector. On rare occasions there can be 
cosmic ray induced events which are unaccompanied by coincident hits in the active 
shield; these arise either from shield inefficiency or from instances wherein a neutron 
or photon with no accompanying charged particles emerged from the cavern walls. 
Contained events of cosmic ray origin are hereafter referred to as "rock events"; these 
can be "shield-tagged" or otherwise of the "zero-shield-hit" variety. 

iv) The "gold events" are data events for which the cavern-liner active shield array 
registered no signal during the event's allowed time window. These events are mostly 
interactions of atmospheric neutrinos B^ . 



C. Processing of events 

Each event included in a contained event sample for physics analysis has passed through 
a standard data reduction chain. Care has been taken to ensure that both data and Monte 
Carlo simulation events pass through the same steps in the chain. For an event to be 
recorded, the central detector requires at least seven contiguous pulse "edges" to occur 
within a 72 /is window before a trigger is registered and the calorimeter modules are read 
out. For MC events, the requirements of the hardware trigger are imposed by a trigger 
simulation code. All events, data and Monte Carlo, are then subjected to a containment 
filter code, which rejects events for which part of the event extends outside the fiducial 
volume. Taken together, the trigger and containment filter requirements will typically reduce 
a sample of MC nucleon decay events by 35% to 65%, depending upon the decay mode being 
investigated. 

Events which survive the trigger and containment requirements, whether they be data 
events (gold or rock) or Monte Carlo (neutrino reaction or nucleon decay), are then subjected 
to two successive "scanning passes" carried out by physicists [0. Each scanning pass 



involves three independent scans. In the first pass, multiprong events are found with an 
overall multiple-scan efficiency of 0.98lo;o4 [^. The first scan pass is designed to ensure 
event quality. Here, residual uncontained events that have not been rejected by the software 
filter are eliminated. Noise events, and events that have prongs ending on detector cracks are 
also rejected at this stage. In the second scanning pass, broad topology assignments are made 
(track or shower - with or without a visible recoil proton, single proton and multiprong) and 
an additional level of scrutiny with regard to containment is provided. Multiprong events 
that survive both scanning passes are then subjected to detailed topology assignment. 

Once topologies have been assigned, events are reconstructed using software tools avail- 
able within our interactive graphics package. Tracks are fitted using a routine that fits a 
polynomial curve through all hits selected and returns the range, energy and initial direc- 
tion vector. There is also a shower processing routine which determines the energy and 
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initial direction vector for clusters of hits that are tagged as showers. The collection of 
reconstructed tracks and showers along with mass hypotheses for the tracks is then written 
into a data summary file, from which the derived quantities such as event energy, invariant 
mass, net event momentum, etc., can be determined. As a measure of event resolution, we 
calculate the difference between Monte Carlo truth and the reconstructed final-state energy. 
The resulting /S.E/ E values for fully- visible final states are (19 ± 2) % for multiprongs with 
prompt muons, and (28 ± 2) % for multiprongs with prompt electrons. These values char- 
acterize reconstruction for nucleon decay (£"*" + hadrons) events and for vi charged-current 
multiprongs near threshold. 

The fractions of event samples which remain after successive application of each of the 
above-mentioned selections, averaged over the full volume of the detector, are given in Table 
II for each nucleon decay mode of this study. For each nucleon decay mode with a particular 
daughter process, a sample of 493 MC events was generated. For the K^ modes, the portions 
of the MC samples which pass all cuts except the kinematic ones are given by the top row, 
leftmost entries of Tables III and IV. 

D. Cosmic ray induced background in "gold" data 

The gold data events analyzed in this work are multiprong events which are mostly 
neutrino-induced, however a small number may arise from rock events which are unaccom- 
panied by hits in the veto shield array. The latter events are initiated by neutrons emerging 
from the cavern rock which are incident upon the calorimeter. In contrast to neutrino inter- 
actions (or nucleon decay) which distribute uniformly throughout the calorimeter volume, 
these events tend to occur at relatively shallow depth into the detector. 

Our estimate of the number of zero-shield-hit rock events in multiprong data is based 
upon a multivariate discriminant analysis of the gold data, neutrino MC, and shield-tagged 



rock samples |3^. For this analysis, discriminant functions characterizing each sample are 



constructed using five event variables (of which several are highly correlated). These include: 
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i) Penetration depth, measured along the event net momentum, ii) vertex distance to the 
closest exterior surface, in) event visible energy, iv) zenith angle of the net momentum 
vector, and v) "inwardness" , defined as the cosine between event pnet and the unit normal, 
pointing inward, of the closest exterior surface of the calorimeter. By fitting the discriminant 
variable distribution from gold data to a combination of neutrino MC and rock discriminant 
distributions, the zero-shield-hit rock contribution is ascertained. Out of the total number 
of multiprong rock events (with and without shield hits), the fraction of events having shield 
hits is / = 0.94 ±0.04. Among 144 gold data multiprongs, 16± 11 events may be cosmic-ray 
induced. Approximately 30% of the zero-shield-hit rock rate can be ascribed to inefficiency 
of the shield array. 

For each nucleon decay channel, we subject the shield-tagged rock multiprongs to the 
same selection cuts which are applied to data multiprongs. The number of zero-shield-hit 
rock events which are background for the channel is then calculated as the product (1 — /)// 
times the number of shield-tagged rock events which pass the cuts. 

E. Neutrino background rates and atmospheric f^ oscillations 

In calculating event rates for neutrino-induced background in a nucleon decay search, 
a consideration arises with the well-known discrepancy between the observed versus pre- 
dicted flavor ratio z/^/z/g for atmospheric neutrinos. In single track and single shower events, 
Soudan 2 observes the ratio-of-ratios {v^/ve observed/expected) to be 0.64 ± 0.11 (stat.) ± 



0.06 (syst.) [^. From the exposure analyzed for this work, we estimate 128 ± 14 contained 
multiprong events to be neutrino events, where the uncertainty is the quadrature sum of 
statistical error and an error from the rock background subtraction. This can be compared 
to the number of neutrino-induced multiprong events predicted by the null oscillation Monte 
Carlo: 154 ±6 events, where the uncertainty represents statistical error from the atmospheric 
neutrino simulation. The difference may be accounted for by atmospheric neutrino oscilla- 
tions u^ to Ux- In any case, evidence for depletion of the atmospheric muon-neutrino flux is 
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sufficiently extensive ^^ that an accounting of the effect in nucleon decay background rates 
is warranted. 

The basis for our atmospheric neutrino background estimates is our reahstic MC simu- 
lation which uses null oscillation fluxes. The disappearance of z/^-flavor neutrinos resulting 
from i/^ -^ Ux oscillations, affects backgrounds initiated by (z/^ + V^) charged-current reac- 
tions; in fact, it reduces this background component. In implementing our correction we 
further assume, as indicated by recent data, that u^ is an active neutrino which is not z/g 
(i.e. Ux = Ur) ^^. Then, the predominant atmospheric neutrino oscillation does not affect 



background arising from charged-current reactions of Ug-fiavoT neutrinos, nor does it affect 
background from neutral currents. Consequently, to correct for z/^-flavor disappearance, the 
number of (z/^ + z7^) charged-current background events estimated from the atmospheric- 
neutrino MC for each nucleon decay channel, has been scaled by the flavor ratio 0.64 ± 
0.13. 

As described below, the neutrino background is very low in the proton decay channels p 
— > /i^K^, K^ -^ TT+TT^ and p -^ /i^K|'. In most other channels, e.g. p -^ yW^K^, K° -^ tt^tt^, 
p -^ e"'"K°, p — > e'^K'i, and n -^ '^K^, K° -^ vr^vr^, the neutrino background is either low or 
otherwise dominated by z/g flavor reactions. In fact, the only decay channel in this work for 
which z/^-flavor disappearance has a significant effect on the lifetime lower limit is neutron 
decay n -^ '^K^, K° -^ tt+tt^. For this particular channel, backgrounds and corresponding 
lifetime lower limits are presented both with and without correction for z/^ oscillations. 

III. CHARACTERIZATION OF NUCLEON DECAY REACTIONS 

For each nucleon decay mode, a Monte Carlo sample is created and processed as de- 
scribed in Section ^, similarly handled are the gold data, the contained shield-tagged rock 
events, and the atmospheric neutrino MC events. These samples are used to determine the 
topological and kinematic properties that differentiate each nucleon decay mode from other 
modes and from the atmospheric neutrino and the rock event backgrounds. 
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A. Event generation 

For each nucleon decay mode, five hundred events are typically generated throughout 
the full detector volume. They are then overlaid onto pulser trigger events which sample the 
detector background noise (from natural radioactivity and cosmic rays) at regular intervals 
throughout the exposure. The effects of Fermi motion within the iron and other nuclei 
which comprise the detector are taken into account in simulation using the parametrization of 
Ref. |^|. For two-body nucleon decays in which the final state momenta would otherwise be 



unique, Fermi motion smears the momenta and thereby complicates final state identification. 
For final states from neutrino interactions and also from nucleon decay wherein pions 
are directly produced, intranuclear rescattering is treated in our Monte Carlo using a phe- 
nomenological cascade model |3^. Parameters of the model, which scale with A, are set 
by requiring that threshold pion production observed in i/^-deuteron {A = 2) and t'^-neon 
{A = 20) reactions be reproduced [^]. For K'' mesons however, inelastic intranuclear rescat- 
tering is expected to be small due to the absence of low-lying K°N {S = +1) resonances 
and is not simulated. 

B. Final state kinematics 

We examined each simulated nucleon decay mode for kinematic characteristics which are 
suitable for selection of data candidates and for separation of backgrounds. Two quanti- 
ties that are generally useful are the invariant mass, Mj„^, and the magnitude of the net 
three-momentum, \pnet\, of reconstructed final state particles. For the distribution in either 
variable, there is an expectation which can be simulated using the nucleon decay MC events. 
For modes in which all final state particles are visible, the invariant mass should fall in a 
range about the nucleon mass (less 8 MeV binding energy) whose width is determined by 
the resolution of the detector. 

If some final state particles are not observable, useful information can still be extracted 
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from the invariant mass. For example, in the case of n ^ z/K^ the mass will correspond to 
that of the reconstructed K^. The nominal 338 MeV/c K*^ momentum from this two-body 
decay will be boosted by the Fermi motion of the decaying neutron and smeared further by 
detector resolution. 

C. Selection contour in the Minv versus |priet| plane 

A scatter plot of invariant mass versus net event momentum can be created for the re- 
constructed final states for each simulation. A region in this plane can be chosen whose 
boundary defines a kinematical selection which can be applied to the data and to the back- 
ground samples. For our searches we consider the application of such a cut to be the 
"primary constraint" and we have formulated a quantitative parameterization of this two- 
variable constraint. We observe that the final-state invariant mass exhibits a distribution 
which is Gaussian to a good approximation. Similarly, the final-state net momentum mag- 
nitude which results from the folding of Fermi motion with the experimental resolution, 
distributes with a shape which is very nearly Gaussian. Consequently, the density distribu- 
tion of points on the invariant mass versus momentum plane can be fitted by a bi-variate 
Gaussian probability distribution function. In principle, when the correlation coefficient is 
set to zero, the form of the fitting function should be 



Zap {x, y) 



1 ~h 



-e 






(3) 



2'K^aaCip 

where o"o and ap are the widths and /Iq and ^p are the means of the two distributions. In 
practice however, the conic section may have semimajor and semiminor axes that are not 
parallel to the ordinate and abscissa of the graph. To allow such a rotation on the plane, a 
new basis (which is a linear combination of the old variables) is constructed with the new 
mean and widths defined on that basis. The form of the fitting function then becomes 



Zap {x, y) = Ae 



'ax + by-fia\^ i fay-bx-fj/}^ ^ 



(4) 
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Note that the two arguments of the exponentials are formulated to be mutually perpendicular 
lines which, when fit, will be coUinear with the major and minor axes of the elliptical contours 
of the sections of the Gaussian surface. 

For each mode in which this representation of a density distribution is used, a lego plot 
for the simulation is created (see Fig. ^). The seven-parameter function above is fitted 
to the simulation using a log-likelihood x^ minimization routine available in the MINUIT 
software package (see Fig. ^d and Fig. |). The bi-variate Gaussian form (Eq. §) was found 
to provide a good fit to the nucleon decay event density distribution in every case, with x^ 
per degree of freedom (dof) values in the range 0.7 to 1.3. 

The projection of the bi-variate Gaussian surface of Fig. 2b onto the Mj„j, versus \pnet\ 
plane is shown in Fig. ^ Here, the shape of the event distribution over the plane is depicted 
using five nested, elliptical contour boundaries. Proceeding outward from the innermost 
contour, the bounded regions contain respectively 10%, 30%, 50%, 70%, and 90% of the 
event sample. From the five regions delineated, we choose the 90%-of-sample contour - 
the outermost, solid curve ellipse in Fig. |^ - for quantitative use in nucleon decay searches 
of this study. In particular, for p -^ ^"""K^ searches we take this outermost contour to 
define our "primary" kinematic selection, namely: A candidate proton decay event must 
have reconstructed (Mj„j,, |p„et|) values which lie within the contour boundary. Concerning 
the interior contours (dashed ellipses of Fig. |^), no quantitative utilization will be made. 
However, by superposing event coordinates onto a contour set such as shown in Fig. |^, a 
degree of insight is afforded as to whether an event sample as a whole exhibits the kinematics 
of nucleon decay. For this reason we display interior contours along with the 90% contour 
in various kinematic diplots of this study. 

D. Other kinematic selection contours 

The kinematic requirement on final-state invariant mass versus net momentum combina- 
tions which the outermost contour of Fig. represents, is a constraint which nucleon decay 
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will satisfy but which background processes need not respect. The background rejection 
thus afforded can be augmented using contours constructed with other kinematic variables. 

From trials with Monte Carlo event samples we find that, with procedures to be described 
below, the final-state lepton prong can be identified with near certainty (in final states 
e+(7r+7r-) and /i+(7r°7r°)) or with useful efficiency (84% for /i+(7r+7r-) and 82% for e+(7r°7r°)). 
With this identification in place, a distinction can be made between the lepton and the K° 
system. (Unfortunately, the neutral track gap produced by the moving K^ is usually not 
large enough to be resolved in the Soudan 2 calorimeter.) Then, the kinematics allowed 
for the Kg subsystem of proton decay can usefully be specified by a contour boundary in 
the plane of K° invariant mass versus K^ momentum. A further kinematic restriction on 
momentum-sharing between the lepton and the K^ can also be formulated (albeit mostly 
redundant with constraints already mentioned) in terms of an allowed contour region in the 
lepton momentum versus K^ momentum plane. 

One might expect that reconstruction of the charged pion tracks in p — > ^"^K^ final state 
with K^ — > TT+TT^ in Monte Carlo samples would on average yield the K° invariant mass. In 
the Soudan detector, however, scattering processes tend to reduce the range of pion tracks. 
As a result, our reconstruction of MC event images which is based upon track range, yields an 
invariant mass distribution with a mean of 461 MeV/c^. Also, the momentum distribution 
of the reconstructed K° is shifted downward by about 50 MeV/c. 

In our i~^K^g and z/K^ searches described below, supplementary requirements are invoked 
using the above-mentioned additional search contours. Their construction proceeds sim- 
ilarly as with the primary constraint contour of Sect. [1IIC|; each contour represents the 



projection of a bivariate Gaussian fit to the distribution of nucleon decay MC events in two 
kinematic variables. A data event, in order to qualify as a nucleon decay candidate, must 
have kinematic coordinates occuring within each contour projection, the boundary of which, 
by construction, includes the coordinates of 90% of the MC nucleon decay sample. 
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IV. NUCLEON DECAY SEARCH 

A. Search for p — »• /^"'"K^ and p — ^ e+K^ 

We now consider searches for two-body decay modes ^"""K^, for which the entire final 
state is observable in Soudan 2. Although the standard SUSY SU(5) GUT models usually 
predict lifetimes for i'^K'^ modes to be two orders of magnitude longer than those for I/K 
modes |^T|, there are other GUT models in which a charged lepton mode could be enhanced 



100 • 



With ^"""K^ modes, the final state will be reconstructed with invariant mass approximating 
the nucleon mass and with event net momentum within a range compatible with Fermi 
motion in iron folded with the detector resolution. The kinematically allowed region in 
the plane of final-state invariant mass versus event net momentum, which is obtained from 
Soudan reconstructions of Monte Carlo proton decay events, is rather similar for all modes 
of this kind. As shown in Figure § and also in Fig. ^, we have used simulations involving 
four separate proton decay modes including K°, to define a "primary" kinematic selection 
contour in the M^v versus \pnet\ plane. That is, for the proton decay searches involving p 
—>■ /i^K° and p —>■ e'^K^, we require all candidate events to have kinematics compatible with 
the outermost contour depicted in Fig. |^a. As indicated by Figs. ||b, |c and ^, the primary 
contour alone provides a very restrictive kinematic selection. 



/^ 



+K0 
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For this mode we examine two topologies: i) three tracks (from K^ -^ tt+tt^) and ii) one 
track plus 3-4 showers (from K° -^ tt^tt^). In the first case, the track which traverses the 
most calorimeter material and which does not have a visible scatter is taken to be the /i"*"; 
the other two tracks are tagged as pions. An illustrative Monte Carlo event of p — ;► /i^K^ 
with the three-track topology, is shown in Fig. ^. For K° —>■ tt^tt" the lone track is assumed 
to be the muon, and the showers are reconstructed as the 7r°7r° system from the K°. To 
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make optimal use of constraining variables, it is useful and also sufficient to invoke here two 
distinct kinematic contours. In addition to requiring the primary contour constraint on the 
final state M^v versus \'Pnet\i we require the "K°" invariant mass versus momentum to fall 
within the contour of Fig. ^ The distributions of each of our four event samples (simulation, 
neutrino MC, rock and data) with respect to the two different kinematic selection contours, 
for both final states corresponding to p ^ /^^K°, are summarized in Figures § and p. In 
Fig. 1^, events that failed the preceeding contour cut are shown using open symbols. A 
candidate event is required to fall within both contours. (Note that events in Fig. ^ can 
populate overflow regions of Fig. |^, e.g. the track-plus-showers event in Fig. ^ is off-scale 
in Fig. |d.) 

As a possible addition or alternative to the contour of Fig. p, event kinematics can be 
evaluated using a contour in the plane of lepton momentum versus K^ momentum. For our 
/i"''Kg search, we flnd no improvement to be afforded. (For e^K° however, such a contour 
provides extra background suppression and we utilize it in that search.) 

For the K^ — ^ tt+tt^ case, the contour cuts are 88% efficient whereas for K° -^ tt^tt^ the 
combined contour cut efficiency is 80%. The application of these kinematic cuts, together 
with the branching ratio, software and scanning efficiencies, results in an overall proton- 
decay survival fraction of (16 ± 3)% for the K° -^ ti^ti^ case and (6 ± 1)% for K^ -^ tt^tt^. 
No event is observed to satisfy all kinematic constraints for either of the vrvr flnal states. 

For the ^~^K.^ flnal state with K° -^ ti^ti^, the estimated neutrino background of 0.6 
events is mostly due to z/g charged- current vr+vr^ production. No neutrino MC events passed 
the cuts for the /i'''K^ flnal state with K° -^ tt+tt"; neutrino-induced backgrounds are 



expected to be low for this case [36 



From our observations of the K^ —>■ tt+tt and ti^tt^ decay modes, we set a proton decay 



lifetime limit using the formula p3|,p^,E2 



/il + /i2 

Here Np(^n) = 2.87(3.15) x 10'^^ protons (neutrons) in a kiloton of the Soudan 2 detector, 
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Tf = 5.52 kiloton years is the full detector exposure (4.41 kty is the fiducial exposure), and 
ej X Bi{K) are the selection efficiencies given in Table II. The /ij are the constrained 90% 
CL upper limits on the numbers of observed events, and are found by solving the equation 

^ . ^ _ ES S:r=o PJ^i^ b, + fi,)P{n,, h + 112) ... 

with the constraint 

A^l /^2 Ej=l fJ'i 

where P{n,fi) is the Poisson function, e~^fi"'/n\, and the bi are the estimated backgrounds. 
With no candidates found for either fi~^K^ mode, we obtain yUi = yU2 = 2.31. The combined 
lower lifetime limit, at 90% CL and with background subtraction, is then t/B > 150 x 

]^Q30 yga2;g_ 

Our limit for p — > fi~^K^ can be compared to the channel limit previously obtained with 
the Frejus planar iron calorimeter: t/B > 64 x 10'^'^ years at 90% CL [^]. The water 
Cherenkov experiments Kamiokande and IMB-3 have reported limits for p -^ yU+K° which 
are more stringent. Lifetime limits for fi~^K^ are summarized below, following the discussion 
of p — > /x^K° searches. 

p ^ e+K" 

We examine events which have the following topologies: i) one shower plus two tracks 
(from K° -^ TT^TT^), or ii) 4 to 6 separate showers (from K° -^ tt^tt*^). For p —>■ e^K^^, 
three kinematic constraint contours are generated for each of the two daughter modes, and 
the data events and both backgrounds are subjected to the cuts thus specified (see sections 
111C| and [111 Dp . Figure ^ shows the primary contour (Mj„^ vs. Ip^netl) for this case; Fig. |^ 



shows the final contour (e^ momentum versus K° momentum) after the K° mass versus K° 
momentum contour constraint has been applied (the latter contour, not shown, introduces 
constraints as in Fig ^.). In case ii) with 4-6 showers in the final state, the prompt e"*" 



shower must be distinguished from the 7 showers originating with decaying vr^'s. In scanning 
of candidate events, the shower prong most hkely to be the positron is chosen based upon 
the topology at the primary vertex. The momentum of the shower chosen by the physicist 
scanner is subsequently compared to the highest momentum among the other showers in 
the event. If the physicists' choice is within 50 MeV/c of being the overall shower of highest 
momentum, it is accepted as the e"^, otherwise the shower of highest momentum is chosen. 
In simulation, this algorithm selects the true e"*' shower in 82% of the events. Again, events 
that have failed a previous cut are plotted using open symbols. Events considered to be 
candidates are required to be inside of all three contours. 

In the case of the e'^K^ final state with K^ -^ vr+vr", the efficiency of all three kinematic 
selections combined is 86%. The overall efficiency times branching ratio is then 15%. The 
number of candidates observed is one event whereas the total expected background is 0.74 
events. Roughly two-thirds of the expected background comes from z/g charged- current single 
charged-pion production events in which the proton is misidentified as another pion. 

Fig. 1^ shows our single e^K^ candidate in the most orthogonal projection which for this 
event is the cathode (Y) versus drift time (Z) projection. A "prompt" electron shower and 
two tracks are seen to emerge from the primary vertex (denoted MARK 1). The mean pulse 
height for either track is too low for a proton assignment but is typical for a pion (or muon) . 
The longer pion track traverses an insensitive module-to-module boundary region, giving rise 
to the apparent gap in the track. The ionization pulse heights on this pion track indicate that 
the vertex at higher F is a secondary scatter, not the primary vertex. These pulse heights 
are smaller (larger) before (after) the proposed scatter point. The kinematics of this event 
is, however, somewhat atypical for proton decay. In the plane of e"*" versus K*^ momentum, 
this event (depicted by the solid triangle in Fig. |^) lies within our search contour, but 
outside the contour which contains 80% of proton decays surviving our selections for this 
mode. It is of interest to consider the "maximum" number of proton decays which could 
appear as in Fig. ^ and which are allowed in our exposure, assuming the limits published 
by Kamiokande, Frejus, and IMB-3 ||26|j2^j2^ to be correct. Using the range of limit values 
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reported for the e"'"K° mode (see below) and the detection efficiency for this channel in the 
Soudan 2 detector (Table II, third row), of order 0.8 to 3.2 occurrences of proton decay 
resembling the Fig. |^ candidate would be compatible with previous limits. 

For the e"'"K° final state with K^ -^ tt^tt^, 75% of all simulated proton decay events 
survive the three kinematic cuts and the overall efficiency times branching ratio for this 
mode is (8 ± 1)%. For this case no data events pass the cuts and 0.63 background events 
are expected. 

Combining the observations from the two separate decay sequences, we obtain for p — >■ 
e+K° a lifetime lower limit at 90% CL of t/B > 120 x 10^° years. A value of t/B > 76 x 10^° 
years was reported for this mode previously by Frejus p^. Limits for e^K° are summarized 
below, following the discussion of p ^ e^Kj' searches. 

B. Search for p -^ fi+Kf and p -^ e+Kf 

For proton decay p — > i~^K^, the K^ daughter decays account for half of the possible final 
states; we now consider p — > £'''K°. In these decays the K° will interact hadronically in the 
Soudan 2 medium before it has time to undergo 37r decay. The event topologies produced 
by the products of the K^ interaction are quite varied but are usually easily discernible. 
The basic topology searched for is that of a single isolated lepton (either /i"*" or e"*"), with 
the appropriate two-body decay momentum, separated from a secondary vertex produced 
by the hadronic products of the K°-nucleus interaction. Additional selections based on the 
"co- linearity" of the event and the characteristic energy visible from the K^ interaction prod- 
ucts are employed to provide discrimination from atmospheric neutrino induced background 
events and also from inelastic interactions of neutrons originating from cosmic ray collisions 
in the cavern rock. 
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We select events whose topology consists of an isolated niuon track which appears sep- 
arated from a hadronic interaction vertex. A Monte Carlo event for this process is shown 
in Fig. ^ This topology selection is applied to each of the four samples of interest, namely 
the p -^ yU^Kl' full detector simulation, the neutrino reactions of the atmospheric u Monte 
Carlo, the shield-tagged rock events, and gold data events. The sample populations thus 
obtained are given in the top row of Table III; the neutrino background rate (column 3) has 
been corrected for z/^-flavor depletion by z/^ -^ u^ oscillations as described in Section |llii^ . 

Kinematic cuts are then applied sequentially; their nature and order are summarized in 
Table III, and the corresponding event populations which survive each cut are listed for each 
of the four samples. We describe below each of the kinematic cuts individually, in the order 
in which they are applied. 

The isolated lepton is imaged as a single track with ionization compatible with that of a 
muon mass assignment. The lone muon track is required to have a length between 35 and 
120 cm. We then require the gap between the event vertex, identified by the start of the fi~^ 
track, and the K^ interaction to be greater than 15 cm. The mean path gap for reconstructed 
proton-decay MC events is 30 cm, however the distribution is broad and extends beyond 
100 cm. Since the secondary interactions exhibit a wide range of topologies with no single 
topology occuring very frequently, we estimate the visible energy of the hadronic final states 
in a crude way by reconstructing entire final states as single showers. That is, the visible 
energy is taken to be proportional to the total number of gas tube crossings ( "hits" ) among 
the associated tracks and showers, with each tube hit requiring about 15 MeV of energy. This 
procedure is carried out using the Soudan 2 shower processor; the processor fits a direction 
vector to the interaction products as well as assigning an interaction energy. In the events of 
the proton decay simulation, this estimated Kj* interaction energy is observed to fall within 
the range 200 to 1400 MeV. This interval therefore represents our third kinematic selection. 

It is to be expected that the K° interaction products will generally travel in the original 
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direction of the K^. We then expect that the vector from the event vertex to the K° 
interaction and the direction vector returned by the shower processor from the K° interaction 
products, point into the same hemisphere. Finally, since p — > fi~^K^ is a two-body decay, the 
fi^ and the Kj' emerge back-to-back in the rest frame of the parent proton. In the laboratory 
frame, Fermi motion can alter this back-to-back configuration. However, in almost every 
event of the simulation, the reconstructed /i"*" direction and the Kj* path gap vector point 
in opposite hemispheres. These angular cuts form the requirement that the event exhibit a 
"co-linearity" such that the /i^, the K°, and its interaction products are all roughly aligned. 



The angular distributions for the proton decay simulation are depicted in Figures |10| and |Tl 
The survival efficiency for proton decay to satisfy the kinematic cuts alone is 53%. Com- 
bining this efficiency with survival rates from the hardware trigger simulation, the fiducial 
volume cut, plus scanning and topology cuts we obtain an overall detection efficiency for p 
— > /i^K° of (12 ± 1)%. Background from atmospheric neutrinos arises mostly from i/^-flavor 
charged- current reactions. With correction for u^ oscillations, a rate of 0.2 background events 
is estimated for the current exposure. A similar background rate is estimated for rock events. 
The 90% CL lower lifetime limit based on the K^ component only is r/B > 83 x 10^° years. 
The p -^ /W^Kj* lifetime limit of 44 x 10^° years was reported previously by Frejus P^ 



From our search results for /i^K^ and /i^K^, we determine a lifetime lower limit for 
proton decay into fi'^K^ (with the K'' ^ K^^ branching fractions included) in the Soudan 
2 experiment: t/B > 120 x 10^*^ years at 90% CL. Limits for /i+K'^ which are similarly 
prepared (90% CL, background subtracted), have been published previously by Kamiokande 
Frejus ||2|], and IMB-3 [|22|; the hmit values are (120, 54, 120) x 10^° years respectively. 



p -^ e+K° 

As in the case for p -^ /i'^Kj', we search for a single isolated lepton prong which appears 
separated from a secondary multiprong vertex. We first require the shower from the prompt 
positron to fall within an energy range of 120 MeV to 500 MeV. The Kj* 15 cm path gap 
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cut is then invoked, followed by the two distinct angle cuts which require an overall co- 
linearity of the shower and the gap and the hadronic interaction secondaries. The numbers 
of background and candidate events which survive these selections are much larger for this 
e'^K^ search than for the fi~^K^ case, as can be seen by comparing the populations tallied 
in columns two through four of Tables III and IV. This situation arises because multiprong 
events, not infrequently, have gamma conversions which are remote from primary vertices 
and thus mimic the topology of p — >■ e"'"K^. 

The interaction energy of the K^ is again constructed by fitting a single shower to the 
interaction products. For e^K°, the selection previously invoked for the K° apparent energy 
in the fi~^K^ search is made more restrictive; here we require 200 to 1100 MeV. The cuts and 
their effects on the simulation, data and background samples are summarized in Table IV. 

For our p -^ e"''K[' search, the efficiency for survival through the kinematic cuts is 50% 
and the e'^K^ detection efficiency is (11 ± 1)%. From our null oscillation atmospheric neu- 
trino sample, 43% of events originate with z/^-fiavor charged- current reactions and receive 
correction for z/^ —>■ u^ oscillations. The total neutrino background (z/g and z/^ flavor, charged 
and neutral currents) is 2.6 events. Our estimate for rock background without accompany- 
ing hits in the active shield is 0.8 events. We find two gold data events which satisfy all 
of our e'^K'i selections. The resulting 90% CL lower lifetime limit for p -^ e~''K[' is then 
t/B > 51 X 10^° years. The 90% CL limit for p -^ e+K° previously reported by Frejus is 



t/B > 44 X 10^° years [|4 



From the above results we determine a limit for proton decay into e"'"K° with subsequent 
K° decay or K° interaction: t/B > 85 x 10^° years at 90% CL. Limits for e"'"K° previously 
obtained by Kamiokande [|6[, Frejus |2|], and IMB-3 |22|, are (150, 60, 31) xlO^" years. 



C. Search for n — > '^K" 



In SUSY GUTs, the mode n -^ z/K" is similar to p — i> z/K"*" in terms of its underlying 
amplitude (see Fig. [^). However, predictions for the relative rates, r(n — » z7^K°)/r(p — 
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z/^K"*") can vary significantly [|,0]- It is possible that amplitude factors conspire to make 
n — i> 77K° predominant, and so a thorough search is warranted. 

For neutron decay into z^K°, K° —>■ tt+tt^, events consisting of two charged (non-proton) 
tracks are expected. For decays leading to K° -^ tt^tt", we expect four (or three) detectable 



showers of modest energies to result from decays of the tt s 



n -^ i/K° K° -^ 77+ 



TV 



For this decay sequence, the kinematic region in the invariant mass versus net momentum 
plane, which contains 90% of the reconstructed neutron decays from simulation, is delineated 
by the outer contour displayed in Fig. |l2]a. Distributions of two-track events (in which 
neither track is a proton) from the atmospheric neutrino MC and from the shield-tagged 
rock event sample are shown in Figs. [TpD and |12|c, where individual events are displayed as 
solid circles. 

The survival efficiency for the uK^ final state with K° -^ tt+tt" through the topology 
and kinematics requirements of Fig. Oa is 53%. Multiplying this by the K^ branching ratio 



(68.6%) and by the survival efficiencies through triggering and containment criteria (Table 
II) we find Soudan's overall detection efficiency for n — > z/K^, K° -^ tt^tt^ to be (17±2)%. In 
Fig. 0d we observe seven events to satisfy the kinematic requirements of Mj^j, and \pnet\ as 
defined by the elliptical contour. Our expectation for background is 6.1 events, of which 5.1 
events are predicted to have been generated by atmospheric neutrinos (without oscillations) 
and 1.0 event from rock. The distribution of data events in Fig. Od is seen to be different 



from one implied by the neutron decay simulation of Fig. |T^a. Out of seven data events 
which occur within the primary contour, only two events fall within a contour which includes 
80% of the nucleon decay sample. 



The mass versus momentum scatter plot of Fig. |T2|b indicates that the neutrino-induced 



background for this nucleon decay mode is significant. Of the reactions for the atmospheric 
neutrino MC events that lie anywhere in Fig. IT2|b, we find 44% to be single-pion or double- 
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pion production events (in which only one of the pions is imaged) generated by z/^ or 17^; 
an additional 42% of these events are misidentified z/^ or V^ quasi-elastics. Electron-flavor 
neutrino reactions with misidentified electrons comprise another 12% of the background. In 
the case that muon neutrinos are oscillating into other-flavor neutrinos our MC prediction 
for the z/^-flavor atmospheric neutrino background is overestimated by approximately 30%. 
Thus, the above 5.1 event prediction for the two-track topologies must be reduced to 3.6 
events giving a revised rock plus atmospheric neutrino total background estimate of 4.6 
events. The lifetime lower limit for this decay sequence with background subtraction is, at 
90% CL, t/B > 40 X 10^0 years. 

n ^ I7K0 K" ^ ttOttO 

In scanning the MC events of neutron decay which proceeds via the daughter process 
K^ -^ Ti^TT^, we find final state four-shower and also three-shower configurations to be the 
most likely topologies. In the following we consider these two cases separately. 

Three-shower final states: Distributions in the invariant mass versus net momentum 
plane for this case are shown in Fig. |13|. Candidate events are required to have kinematics 
corresponding to the occupation of the elliptical domain depicted in Fig. [l3|a. The detection 
efficiency for n — > VK^, K^ -^ tt^tt^ -^ 3 showers is 3.3%. The number of neutrino plus rock 



background events implied by Figs. 13b and 13c is 3.4 events. Interrogation of the atmo- 



spheric neutrino MC events indicates that the neutrino background reactions are generated 
by z/e and z7e multi-pion production charged-current events. For the cases in which charged 
pions were produced, intranuclear charge exchange forced the observed final meson state 
into a 7r° from which the proton decay products were reconstructed along with the prompt 
electron shower. Fig. p^ shows that seven candidate events are actually observed. From 
the published Kamiokande limit for this channel together with the low detection efficiency 
of this experiment (Table II, row 8), we infer that the signal into three showers in Soudan 
2 should be 0.3 events or less. Thus, the possibility that our observed 3.6 event excess 
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could be nucleoli decay is disfavored. We set a lifetime lower limit for this channel which is 
t/ B > 5.6 X 10^° years. 

Four-shower final states: Distributions in Minv versus \pnet\ for four-shower events of 
the neutron decay simulation, of the atmospheric neutrino Monte Carlo, of the rock event 
sample, and of the gold data sample, are shown in Figs. Ilia through [Tld respectively. The 



overall selection efficiency for this case is 4.9%. The number of candidate events observed is 
two (see Fig. [T^) whereas the total expected background is 1.1 events. For this four-shower 
case, the background comprising the atmospheric neutrino Monte Carlo sample is composed 
of Ue and z7e multi-pion production charged- current events in which the prompt electron 
shower and showers from the resulting vr^s make up the four reconstructed prongs. Using 
four-shower events only, the lifetime lower limit at 90% CL is t/B > 18 x 10^° years. 

Given the absence of a significant overall signal pattern, we have combined all three 

limits from the separate cases treated above. For n -^ z7K° we obtain a limit of 51 x 10^° 

years. For n —>■ z7K° then, the lifetime limit at 90% CL from our experiment is 26 x 10^° 

years. Limits for n — > I7K° previously published by Frejus [Q, Kamiokande ||2^, and lMB-3 

2^ are (15, 86, 30) x 10^*^ years respectively. 



V. STATISTICAL AND SYSTEMATIC ERRORS 

In the lifetime lower limit determinations of this work, there are statistical and systematic 
errors which arise in event detection and final state recognition for each decay mode, and 
there are errors, primarily systematic, which arise from background estimations. We review 
here the variations which may be introduced by error sources, and estimate the uncertainty 
^'Tn /'Tn on the lifetime limits obtained. 

Table II displays the survival efficiencies for each of ten distinct nucleon decay plus 
K decay sequences, through the selections imposed by triggering, filtering, scanning and 
kinematic cuts. Knowledge of these efficiencies is limited by the statistics and systematics 
of the simulation used for each sequence. Combined statistical errors are given with the 



26 



total efficiency of eacli sequence in tlie riglit-most column of Table 11. These errors range 
between 11% and 18%; as can be seen from Eq. (|^), they contribute directly to Atj^/t^. 
Systematic error contributions to detection efficiencies could arise through inaccuracies in 
the nucleon decay simulation, including kinematic cut inaccuracies in simulation versus data. 
With regard to the latter, we note that the Monte Carlo has been tested against data taken 
at the Rutherford Laboratory ISIS test beam facility |^3|. The comparisons were carried 



out using beams of electrons, muons, and pions with a variety of momenta extending up to 
400 MeV/c. For particle energies relevant to nucleon decay, the comparisons are reassuring 
and would seem to exclude significant kinematic offsets. 

Concerning nucleon decay simulation, an uncertainty arises from inelastic intranuclear 
rescattering of K° within parent nuclei. A loss of 10% of K° born in oxygen, as a result of 
absorption, charge exchange, and large-angle scattering, has been estimated by 1MB [p!7| . 
Given the uncertainties in such an estimation we prefer not to include an explicit efficiency 
factor, however we conclude that losses of up to 15% in the Soudan medium are possible. 

Errors also enter the lifetime lower limit calculation (Eqs. (|§)-(0)) through the estimates 
of backgrounds which include atmospheric neutrino events and the cosmic ray muon induced 
rock events. Uncertainties in backgrounds have little effect on our limits for p -^ /i'''K° and 
p — > fi^K^ wherein zero candidates are observed and low background is expected. This is 
not the case for n —>■ z/K^ however, for which the numbers of candidates and background 
events are sizable. The errors in our rock background estimates arise primarily from the 
finite statistics of shield-tagged rock control samples. Statistical errors arise from the finite 
size of the atmospheric neutrino Monte Carlo sample which corresponds to an exposure of 
24.0 fiducial kiloton years. For the individual e"'"K^ and z/K° processes (rows 3, 4 and 6-10 
in Table II) the above-mentioned uncertainties imply At^/tn — 15-20%. 

Background rates for neutrino reactions are based upon our atmospheric neutrino Monte 
Carlo, and so the estimates are subject to uncertainties in the absolute neutrino fiuxes 
(20%), in neutrino cross sections (30%), in the treatment of intranuclear rescattering (30%), 
and in the treatment of other nuclear effects, e.g. Fermi motion and Pauli blocking (20%) 
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4^ . Inherent to our correction of (z/^ + z7^) charged- current background rates to account for 
neutrino oscillations, is the uncertainty (20%) in the Soudan 2 flavor ratio. Uncertainties 
in the backgrounds, although sizable, affect the lifetime lower limit calculation through 
variation of the 6j of Eq. (|]); the variation propagates to the /Xj of Eq. (|^) via solution of 
Eq. (^) with constraint (0) - a non-linear relation. Propagating the above listed uncertainties 
through for the n -^ z/K° case, we find Atn/tn — 16%. 

We conclude that the uncertainty Atn/tn on the lifetime lower limits reported here 
arising from the sources described above is approximately 20% for proton decay into /i'''K° 
and fJ^'^K'i, and may be as large as 32% for e"'"K^ and e~^K^ final states, and for neutron 
decay to z/K°. 

VI. DISCUSSION AND SUMMARY 



We report results from a search for nucleon decay into two-body, lepton-plus-K" final 
states. Current SUSY Grand Unified Theories propose that these nucleon decay modes 
occur and may be predominant (Refs. 7-12); our investigation probes the lower range 
of published lifetime predictions. The search is carried out using Soudan 2's fine-grained 
tracking calorimeter of honeycomb lattice geometry, a detector which images non-relativistic 
as well as relativistic charged tracks and provides dE/dx ionization sampling. The fiducial 
exposure utilized, 4.41 kiloton- years, is twice as large as any previously reported using the 
tracking calorimeter technique. For each of the decay channels investigated, selections have 
been designed which reduce the neutrino and cosmic ray background to levels of few events, 
while maintaining sufficient detection efficiency to allow a sensitive search. In each of the 
channels we find zero or small numbers of nucleon decay candidates; the latter occur at or 
below rates calculated for background. No evidence for a nucleon decay signal is observed, 
and we report lifetime lower limits r/ B at 90% CL as summarized in Table V. In the Table, 
the limit r/ B given for each K° channel is larger than the limit for the corresponding K° 
channel due to the branching fraction B{K^ -^ K°) = 0.5. For the K° channels, the listed 
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t/ B values represent our best measured limits obtained by summing the K^ and K^ channels. 
Under the assumption of CP invariance in the decay, i.e. i?(K° — > K°) = i?(K° -^ K°) = 0.5, 
better limits equal to the K^ limits could be inferred. 

For the p -^ /U^K° and p -^ yU^K° channels we estimate the background to be less than 
one event and we find no candidate events. The lifetime lower limits at 90% CL thus obtained 
for (/i'^Kg, /i~'"K[') are (150,83) x 10^° years. Previous limits for these channels, reported 
by Frejus [H and hsted by the Particle Data Group (PDG) g||, are (64,44) x 10^° years. 



By combining these channel limits, the limit for p — > /i^K*^ in the Soudan 2 experiment is 
obtained: r/ B > 120 x 10'^° years. Previous limits for this mode from Kamiokande, Frejus, 
and IMB-3 are (120, 54, 120) x 10^° years respectively [HHH . 



For p -^ e+K° and p — >■ e"'"K5' channels we estimate background rates of one event and 
3.5 events, and we observe one candidate and two candidates respectively. Our lifetime 
lower limits for (e"'"K^, e"'"K°) are (120,51) x 10'^° years, an improvement upon the limits 
(76,44) X 10^° years established previously by Frejus [^,Q. For the mode p — > e"'"K°, we 
obtain r/ B > 85 x 10^° years, to be compared with previous limits by Kamiokande, Frejus, 
and IMB-3: (150,60,31) x 10^° years respectively [|6|,||,|2|. 



For neutron decay into z/K^ we observe candidates appearing at rates which are com- 
patible with background estimates. Our lifetime lower limit for z/K° is t/ B > 26 x 10^° 
years. Although the Soudan 2 limit for this mode is the most restrictive ever obtained using 
the tracking calorimeter technique, limits of (86, 30) x 10^° years have been reported by 
Kamiokande and IMB-3 ||,|22|. 
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TABLES 



Model 



Leading Mode(s) Predicted Lifetimes Authors [Ref.] 



SUSY SU{5) n -^ VK^, p -^ VK+ ~ lO^^y Hisano, Murayama (1993) 



Discrete (S^)^ 



p -^ e+K° 



lO'^^y Carone et al. (1996) H 



SUSY 50(10) with 
"reahstic masses" 



p -^ l+K^ 



Babu & Barr (1996) |l 



SUSY 50(10) n -^ vK^, p -^ i/K+ ~ lO^^y i^^^^ ^ j^^^j (I997) 



Anomalous U{\) 
from super strings 



p -^ /u+RO < 1032y 



Irges et al. (1998) P 



SUSY 50(10) with p ^ ^+K0, ^+7r° 



see-saw v masses 



Babu et al. (1998) |12| 



TABLE I. SUSY model predictions for nucleon decay into 'lepton + K final states. 
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Decay 

Mode 



Daughter Hardware Contain- Event Topology Kinematic BR ex BR 

Process Trigger ment Quality Selection Cuts 

Filter Scans 



p - 


-.^+ko 


K°- 


-^ vr+TT 


0.97 


0.67 


0.72 


0.56 


0.88 


0.69 


0.16±0.03 


p - 


-M+KO 


K°- 


-^^V 


0.99 


0.57 


0.67 


0.63 


0.80 


0.31 


0.06±0.01 


p- 


-.e+KO 


K°- 


-^ 7r+7r~ 


0.97 


0.68 


0.73 


0.51 


0.86 


0.69 


0.15±0.03 


p - 


-.e+K" 


K°- 


-^ttOttO 


0.99 


0.61 


0.70 


0.78 


0.75 


0.31 


0.08±0.01 


p- 


-.^+K0 


K? 


^ int. 


0.99 


0.59 


0.77 


0.53 


0.53 


1.0 


0.12±0.01 


p- 


.e+KO 


K? 


^ int. 


0.97 


0.61 


0.74 


0.48 


0.50 


1.0 


O.llibO.Ol 


n - 


-.I7K0 


K°- 


-^ TT^vr^ 


0.87 


0.73 


0.72 


0.56 


0.94 


0.69 


0.17±0.02 


n - 


-.I7K0 


KO^. 


tOttO (3 7) 


0.93 


0.68 


0.61 


0.31 


0.90 


0.31 


0.03±0.01 


n - 


-.I7K0 


K°-. 


IT^TT^ (4 7) 


0.93 


0.68 


0.61 


0.44 


0.94 


0.31 


0.05±0.01 



TABLE II. Event sample survival fractions (columns 3 through 7) and channel detection effi- 
ciencies (column 9) for nucleon decay final states containing K" mesons. The events simulated for 
these Monte Carlo samples were generated within the full volume of the detector, consequently the 
survival rates reported in the Containment and Event Quality columns include removal of events 
which extend outside the fiducial volume. 
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Event 


p -^ fi+Ki 


Atm. u 


Rock 


Data 


Selections 


Simulation 


Background 


Background 


Events 




(493 evts) 


(4.41 fid.kty) 


(4.41 fid.kty) 




Track-Gap-Interact. Topology 


110 


2.2 


0.8 


1 


35 < l^+ < 120 cm 


105 


2.0 


0.6 


1 


15 cm< KJ* path gap 


80 


0.9 


0.5 


1 


200 < K^ "energy" < 1400 MeV 


69 


0.4 


0.4 





cos(K° • interaction) > 


63 


0.4 


0.3 





cos(K[' • /u+) < 


58 


0.2 


0.2 






TABLE III. For the p -^ ^^KJ* search, the evolution of event sample populations resulting 
from successive application of topology and kinematic selections. The four columns to the right 
of each Event Selection show survival populations for the nucleon decay simulation, for neutrino 
background events of the atmospheric ly Monte Carlo, for the rock background (estimated from 
shield-tagged data), and for gold data events. The atmospheric neutrino and rock background 
populations are normalized to the 4.41 fiducial kiloton-year exposure. 
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Event 


p ^ e+ RO 


Atm. u 


Rock 


Data 


Selections 


Simulation 


Background 


Background 


Events 




(493 evts) 


(4.41 fid.kty) 


(4.41 fid.kty) 




Shower-Gap-Interact. Topology 


101 


17.5 


3.2 


25 


120 < Pe+ < 500 MeV/c 


92 


17.0 


3.1 


25 


15 cm < K^ path gap 


67 


14.0 


2.7 


18 


200 < K^* "energy" < 1100 MeV 


61 


6.7 


1.8 


8 


cos(K° • interaction) > 


55 


2.9 


0.9 


3 


cos(K[' • e+) < 


51 


2.6 


0.8 


2 



TABLE IV. For the p -^ e'''K^ search, the evolution of event sample populations resulting from 
successive application of topology and kinematic selections. Rates for surviving atmospheric neu- 
trino and cosmic-ray neutron-induced backgrounds are for the 4.41 fiducial kiloton-year exposure. 
The net backgrounds (bottom row, columns two and three) are higher than those obtained for 
/i+K^" (see Table III). 
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Decay Mode 



Final State 



exB.R. 



^Bk 



Total^^ Data t/B x IQ-^'Jy 



P ^ /^+kc 



fl' TT ' TT 



/X+TfO^O 



0.16 < 0.2 <0.2 

0.06 0.6 0.6 



150 



e+ K^ 



e^TT^Tr 



e+TrOyrO 



0.15 
0.08 



0.6 
0.4 



0.7 
0.6 



1 




120 



P - Z^+K? 


K? 


-^ interaction 


0.12 


0.2 


0.4 





83 


p ^ e+KO 


K? 


-^ interaction 


0.11 


2.6 


3.5 


2 


51 


p ^ ;U+K0 


fJ- 


+ (KO + KO) 


0.17 


<0.9 


<1.2 





120 


p ^ e+KO 


e" 


^KO + KO) 


0.17 


3.5 


4.9 


3 


85 


n^r^KO 




TT+TT" 


0.17 


3.6(5.1) 


4.6(6.1) 


7 


51(59) 


(three - showers) 




ttOttO 


0.03 


2.6 


3.4 


7 




(four - showers) 




ttOttO 


0.05 


0.6 


1.1 


2 




n^i^RO 


KKO + KO) 


0.13 


6.8(8.3) 


9.1(10.6) 


16 


26(29) 



TABLE V. Background-subtracted hfetime lower limits at 90% confidence level from Soudan 
2. Correction of neutrino background for z^^-flavor depletion by oscillations has significant effect 



for n -^ J^K^; values without this correction are given in parentheses. For K channels, we list our 
best measured limits obtained by summing the K^ and K^ channels. 
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FIG. 1. A nucleon decay amplitude in SUSY GUTs. 
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FIG. 2. a) Simulation for decay mode p — > ^"^K^, displayed as a lego plot in Mj„^ vs. |pnet| 
space, b) The bi-variate Gaussian surface fitted to the distribution in Fig. ga. 
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FIG. 3. The contours of the bi-variate Gaussian surface of Fig. g projected onto in the Minv 
vs. \pnet\ plane. The region bounded by the outermost contour contains 90% of reconstructed 
proton decay events and defines the primary kinematic selection used to identify proton decay 
candidates. 
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FIG. 4. For proton decay modes p — > ^"""K^, the "primary" kinematic selection contour (out- 
ermost contour) together with event distributions, in the Minv versus \pnet\ plane. Distributions 
show a) the proton decay simulations, b) atmospheric neutrino MC events, c) rock events, and d) 
data events. The fJ-'^K.^ (e"^K^) final states are depicted using solid circles and squares (triangles 
and stars). 
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FIG. 5. Monte Carlo event with full detector response for proton decay p -^ ^"^K^, K^ — > vr+vr". 
Here, the anode {X) versus drift time {Z) projection has been selected from the three scanning views 
which display anode-time, cathode-time, and anode-cathode images. The open square denoted 
'MARK 1' depicts the reconstructed primary vertex. 
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FIG. 6. For the K^ subsystem of p ^ ^^K^ candidates, event distributions plus kinematic 
selection contour in the plane of Minv versus Ipk^ I • Distributions show a) proton decay simulations, 
b) atmospheric neutrino MC events, c) rock events, and d) data events. Open symbol events fall 
outside the primary contour of Fig. 0. 
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FIG. 7. The plane of e"*" momentum versus "K^" momentum for p -^ e+K^ candidates, showing 
event distributions plus kinematic selection contours. Distributions show a) proton decay simu- 
lations, b) atmospheric neutrino MC events, c) rock events, and d) data events. Open symbols 
denote events that fall outside of the primary contour and/or the K^ contour of the K^ subsystem 
contour. 
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FIG. 8. Data candidate for proton decay p — > e'''K^, Kg — > 7r"'"7r , shown in the cathode (1") 
versus drift time (Z) projection. 
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FIG. 9. Monte Carlo event for p -^ Ai^KJ' with subsequent Kj' interaction in the calorimeter 
medium, shown in the cathode versus drift time projection. 
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FIG. 10. Angle between the K.^ flight path and visible net momentum of its subsequent 
hadronic interaction. 
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FIG. 11. Angle between the /U"*" track and the Kj' flight path, for p -^ ^"'"KJ'. 
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FIG. 12. For n -^ i^K°, K^ -^ vr+vr^, event distributions and kinematic selection contour in 

the Minv versus \'Pnet\ plane. Distributions show a) the neutron decay simulation, b) atmospheric 

neutrino MC events, c) rock events, and d)the data events. Background events and data candidates 

having the two-track topology of this mode, are shown via solid circles. 
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FIG. 13. For n -^ i^K^, K^ -^ ir^ir^ yielding three visible showers. Distributions in the Minv 
versus \pnet\ plane show three-shower events of a) the neutron decay simulation, b) the atmospheric 
neutrino MC, c) the rock events, and d) the data events. 
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17K3, K*^ -^ TT^TT^ yielding four visible showers. Distributions in the Minv 
versus \pnet\ plane show four-shower events of a) the neutron decay simulation, b) the atmospheric 
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